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A  COHERENT  NONLINEAR  THEORY  OF  AURORAL  KILOMETRIC  RADIATION: 

II.  DYNAMIC  INTERACTIONS 


I.  Introduction 

In  a  recent  paper  [Grabbe,  et  al.  1980] ,  a  theory  of  auroral 
kilometric  radiation  was  proposed,  in  which  electromagnetic  noise  is 
amplified  by  interaction  with  low  frequency  coherent  quasineutral 
density  fluctuations  created  by  electrostatic  ion  cyclotron  (EIC) 
waves,  in  the  presence  of  precipitating  auroral  electron  beams.  The 
result  is  a  three-wave  parametric  process  in  which  a  beat  wave  is 
produced  that  can  interact  with  the  beam,  much  like  the  theory  of 
Palmadesso,  et  al  [ 19 76 J .  It  was  found  that  when  the  wave  frequency 
is  in  the  right  range,  the  electromagnetic  wave  is  negative  energy  in 
the  rotating  frame  of  the  beam  electron  and  undergoes  a  convective 
instability.  The  basic  requirements  for  the  instability  were  found 
to  be: 

(1)  Minimum  beam  density: 


n  \  k  (Av) 2 

_b_  >  z - 

n  /  2w  v, 

o'  ce  b 


(1) 


where  n.  and  a  are  the  beam  and  plasma  density,  respectively,  v,  and 
bo  o 

Av  the  beam  velocity  and  thermal  spread  in  velocity  space,  k^  the 
wave  vector  component  along  the  magnetic  field,  and  the  electron 

cyclotron  frequency. 

(2)  Accessibility  to  free  space  (w  >  coR  where  OJg  is  the  right  hand 
cutoff) : 


oj  2  <  k  v,  iii  (2) 

pe  z  b  ce 
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(3)  Frequency  range: 


oj  <  oj  +  k  v. 
ce  zb 


(3) 


Condition  (1)  typically  requires  nb  >  10-3  nQ  and  condition  (2)  typi¬ 
cally  requires  local  depletion  of  the  plasma  density  such  that 
Wpe  <  ^*2  “>ce*  The  latter  is  in  very  good  agreement  with  the  observed 
density  depletion  [Benson  id  Calvert.  1979]. 

Combining  conditions  (1)  -  (3)  gives  the  following  limits  on 
the  frequency  and  propagation  directions  of  the  wave  for  amplification 
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Eq.  (4)  gives  radiation  in  a  narrow  frequency  just  above  the  right 
band  cutoff,  while  Eq.  (5)  normally  limits  the  propagation  to  be 
almost  (but  not  quite)  perpendicular  to  the  magnetic  field  [Grabbe, 
1980].  Furthermore,  the  0-mode  has  no  such  unstable  frequency  range, 
hence  the  X-mode  is  the  predicted  polarization.  All  of  these  pre¬ 
dictions  are  in  excellent  agreement  with  observation. 

The  above  conclusions  were  based  on  a  steady  state  model,  in 
which  amplitude  of  the  density  fluctuations  was  assumed  to  be  approxi¬ 
mately  constant.  This  is  valid  if  the  energy  in  the  density  fluc¬ 
tuations  is  replenished  by  the  beam  or  other  sources  at  approximately 
the  same  rate  as  it  is  being  used  up.  However,  the  Feynman  diagram 
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for  the  three-wave  process  (Fig.  1)  reveals  that  a  more  dynamical 
process  is  taking  place.  Not  only  is  energy  being  resonantly  trans¬ 
ferred  from  the  density  fluctuation  to  the  electromagnetic  wave  in  the 
appropriate  frequency  band  because  of  the  beam,  but  the  energy  the 
beam  injects  into  the  beat  wave  is  being  transferred  back  to  the 
density  fluctuations  and  the  electromagnetic  wave  because  of  a  finite 
three  wave  coupling  coefficient.  This  coupling  coefficient  was  ignored 
in  the  steady  state  theory,  but  must  be  included  to  understand  the  full 
dynamical  process. 

The  purpose  of  this  paper  is  to  study  the  dynamics  of  the  full 
three  wave  process,  including  temporal  variations  in  the  density  fluc¬ 
tuations  and  a  finite  three-wave  coupling  coefficient,  in  order  to 
confirm  the  predictions  of  the  steady  state  tneory.  In  Sec.  II  we 
introduce  a  set  of  coupled  nonlinear  rate  equations  for  the  evolution 
of  each  of  the  three  waves  in  the  Feynman  diagram,  and  discuss  steady 
state  solutions.  Several  numerical  solutions  of  these  equations  are 
presented  and  discussed  in  Sec.  Ill  for  various  typical  parameters. 

The  principle  conclusions  are  summarized  in  Sec.  IV. 
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Fig.  1  —  Feynman  diagram  for  the  three-wave  interaction  involved  in  the  amplification  of 
radiation  to  produce  auroral  kilometric  radiation.  The  processes  numbered  are  (1)  induced 
absorption,  (2)  beam  amplification,  (3)  induced  (and  spontaneous)  emission,  and  (4)  reso¬ 
nant  amplification.  The  rate  coefficient  for  each  process  is  shown. 
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II.  Dynamical  Model 


The  dynamical  processes  involved  in  the  amplification  process 
which  produces  the  AKR  are  summarized  in  the  three-wave  Feynman  diagram 
in  Fig.  1.  These  processes  are  (1)  Induced  absorption  of  the  X  mode  by 
the  density  fluctuations  to  create  the  quasimode  or  beat  wave  (2)  Beam 
amplification  of  the  beat  wave  (3)  Spontaneous  and  induced  emission  of 
the  X-mode  by  the  beat  mode  (4)  Resonant  interaction  between  the  density 
fluctuations  and  the  X  mode  in  which  energy  is  transferred  to  the  X  mode 
The  conditions  for  the  last  process  were  given  in  the  introduction,  and 
the  resulting  spatial  growth  rate  was  given  by  Eq.  (28)  in  Grabbe,  et  al 
[1980]: 


We  want  to  formulate  the  rate  equations  for  the  aforementioned 
processes  in  terms  of  the  (quantum)  occupation  number  density  of  the 
waves. 


N.  = 
1 


I E  I  ^  /  8ttuj. 
1  J  1 


(7) 


where  t  is  the  wave  electric  field.  We  designate  Nq  as  the  number 
density  of  the  X-mode,  for  the  density  fluctuations,  and  for 
the  beat  wave.  We  then  have  the  following  contributions  to  the  rate 
equation  for  each  process  numbered  in  Fig.  1: 
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(1)  Induced  absorption. 

dN, 


b 

dt 


-  W  N  N. 
o  1 


(8a) 


“dt 


W  N  N, 
o  1 


(8b) 


dN 

— -  =  2W  N  N. 

dt  o  1 


(8c) 


where  W  is  the  three  wave  coupling  coefficient  [Tsytovich,  1974] 


W  = 


2  3  2 

e  co  .  to  m . 
ei  pe  l 


o  3  3  4,2 

87rme  co  vTe  k 
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(2)  Beam  amplification. 

dN, 

— —  =  r  n 

dt  ‘bb 


(10) 


where  we  take  the  beam  driven  growth  rate  of  the  quasimode  to  be  the 
usual  form  for  a  beam  plasma  instability  [Briggs,  1964]: 

,1/3 
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(3)  Spontaneous  and  induced  emission. 

dN 

-7 7  =  W  N.  (1+N  +N  ) 

dt  b  o  1 


(12a) 
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-rf  =  WN.d+N+N.)  (12b) 

dt  b  o  1 

dN 

— —  =  -2W  N,  ( 1+N  +N  )  (12c) 

dt  b  o  I 

The  spontaneous  emission  term  (first  term  on  the  right  hand)  is  normally 
negligible,  since  normally  Nq  >>  1  or  >>  1  (N  and  N^ ,  taken  as 
dimensionless  occupation  numbers). 

(4)  Resonant  interaction. 


dN  . 

— r2-  =  yN  N1/ 

dt  '  o  1 


(13a) 


=  -yN  N^/2 
dt  '  o  1 


(13b) 

Here  the  dependence  of  the  resonant  growth  rate  f  of  the  X-mode  on  the 


density  fluctuation  N^  has  been  explicitly  factored  out: 


YN1/2 


(14) 


Here  the  temporal  growth  rate  can  be  expressed  in  terms  of  the  spatial 
growth  rate  by  multiplying  by  the  group  velocity 

r  =  2  k  c  (15) 

where  the  factor  2  represents  the  conversion  between  electric  field 
growth  rate  an  the  quantum  density  growth  rate.  Thus  from  Eq .  (6) 


at  its  maximum  value. 

Combining,  (1),  (2),  (3),  (4)  we  find  the  complete  set  of 
equations  to  be 


(16) 
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(17a) 


Since  the  steady  state  model  only  considers  the  contribution  of 
process  (A),  it  only  gives  one-half  of  the  growth  rate. 
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III. 


Numerical  Results 


To  obtain  typical  values  for  the  growth  rates  and  ooupling 

coefficient  we  use  Eqs.  (9),  (11)  and  (16).  If  we  take  typical  source 

n,  -3  n, 

region  values  n,>10  n  ,  u)  <0.2  u>  ,  n,n.O. 3-0.5  n  ,  we  find  the  following 
&  b  o  pe  ce’  1  o’ 

normalized  values  when  the  N  's  are  normalized  to  dimensionless  values 
-  a 

T,  >  6  x  104  sec  *  (19a) 

D 

T  =  y  N^>  6  x  104  sec  *  (19b) 

W  n,  0.05  sec  ^  (19c) 

Here  N  =  (e^/m  ^c^u)N  ,  y  =  (e/m  c/w  )  and  W  =  (m^c^/e^)W.  These  values 
a  e  a  e 

will  be  used  as  a  guide  for  input  to  the  numerical  calculations. 

The  kinetic  equations  in  Eq .  (17)  were  solved  for  several  values 
of  the  growth  rate  and  initial  conditions,  although  the  initial  value 
of  was  always  set  to  the  normalized  value  on  1 .  A  classical  Runge-Kutta 
integration  routine  was  used  initially  but  proved  to  be  too  inefficient, 
so  it  was  replaced  by  a  stiff  integration  scheme  called  CHEMEQ  [Young, 
1980],  A  sampling  of  the  results  is  shown  in  Figs.  2-4.  Included  in  the 
graphs  is  the  ratio  of  the  time  dependent  growth  rate  F*(t)  to  the 
resonant  growth  rate  T,  where 
dN 

—£■  =  r*(t)N  .  (20) 

at  o 

The  graphs  show  that  there  are  two  principle  stages  of  growth 
of  the  X-mode .  In  the  initial  stage  the  growth  rate  is  just  the 
resonant  growth  rate  dNo/dt=yNoN^ .  The  reason  for  this  stage  is  that 
not  much  energy  has  been  transferred  into  the  beat  wave,  so  all  of  the 
X-mode  energy  is  coming  from  the  density  fluctuations.  This  stage 
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-  y* 

(105  sec-1) 


Fig.  4  —  Graph  of  the  ratio  of  the  growth  rate  of  the  wave  as  function  of  time  to  that 
found  in  the  steady  state  theory,  for  the  conditions  in  Fig.  3.  Note  that  it  is  normally 
2n,  for  various  non-negative  integer  values  of  n. 


normally  saturates  when  most  of  the  energy  in  the  density  fluctuations 


is  absorbed.  The  second  stage  starts  when  a  significant  amount  of  energy 
has  been  pumped  into  the  beat  wave  by  the  beam,  and  this  energy  begins 


being  transferred  to  the  density  fluctuations  and  the  X-mode .  This 

* 


second  stage  grows  at  dNo/dt=2YNQN£,  or  twice  the  resonant  growth  rate. 
The  graphs  of  f*(t) /I'=Y*(t) /y  shows  thc’-e  may  also  exist  subsequent 


1,  W 

stages  in  which  the  growth  rate  occurs  at  dNQ/dt=4YNoN  “ ,  dNo/dt=8YNQN^ 


etc.  Then  the  process  reaches  a  final  level  of  saturation  when  almost 
all  of  the  energy  of  the  density  fluctuations  and  the  beat  wave  has  been 


depleted.  The  stages  of  growth  confirm  that  the  rate  of  growth  is 


governed  by  the  resonant  growth  rate  T=yN^,  as  predicted  by  the  steady 
state  theory. 

A  comparison  of  the  saturation  amplitude  of  the  X-mode  shows 
that  it  always  is  several  orders  of  magnitude  above  the  initial  value 


of  the  density  fluctuations  N^(t=0) .  This  shows  that  virtually  all  of 


the  energy  comes  from  the  beam,  rather  than  from  the  initial  level  of 
the  density  fluctuations.  Furthermore,  when  is  varied  on  sample  runs 
with  all  other  parameters  being  held  constant,  then  the  saturation 
amplitude  increases  as  F^  is  increased;  thus,  F^  determines  the  total 
energy  the  electromagnetic  wave  can  absorb.  A  statistical  analysis  of 


the  saturation  level  N  _  from  many  sample  runs  reveals  the 

o  sat 


scaling  law 


n  .  n,  rL 

o  sat  b 


4/3 


(21) 


From  Eq.  (11),  we  have  the  dependence  on  the  relative  beam  density 
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A  derivation  of  this  scaling  law  is  given  in  the  Appendix. 

A  comparison  has  been  made  on  sample  runs  for  different  initial 
values  of  all  other  parameters  held  constant.  It  is  found  thit 
although  the  growth  rate  of  the  X  mode  increases  linearly  with  N^,  the 
saturation  amplitude  is  relatively  independent  of  the  initial  value  of 
N^!  The  significance  of  this  for  the  growth  process  in  the  finite 
density  depleted  cavity  introduced  in  the  steady  state  model  [Grabbe , 
et  al.,  1980]  is  the  following.  If  the  X-mode  saturates  before  propagating 
out  of  the  growth  region  in  the  cavity,  then  an  increase  in  the  initial 
levels  of  the  density  fluctuations  does  not  have  a  very  important  effect 
on  level  of  AKR  produced.  However,  if  the  X-mode  propagates  out  of  the 
growth  region  before  reaching  its  saturation  level,  increasing  the 
Initial  level  of  the  EIC  density  fluctuations  would  bring  the  AKR  closer 
to  its  saturation  value.  Calculation  done  in  the  steady  state  model 
would  suggest  the  former  case  occurs  more  often  than  the  latter. 
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IV.  Summary  and  Conclusions 

We  have  formulated  the  dynamics  of  the  three-wave  process 
involved  in  amplification  of  the  X-mode  to  produce  AKR  in  terms  of  a 
set  of  coupled  rate  equations.  An  analysis  of  these  coupled  rate 
equations  has  confirmed  the  conclusions  drawn  from  the  steady  state 
theory.  It  was  found  that  the  growth  rate  was  determined  by  the 
resonant  growch  rate  found  in  the  steady  state.  The  growth  was  seen 
to  occur  in  stages:  an  initial  stage  in  which  the  growth  rate  was 
just  the  resonant  growth  rate,  and  subsequent  stage  of  growth  at  2n 

times  the  resonant  growth  rate,  where  n=l,2,3 .  This  is  followed  by 

a  saturation  of  the  wave. 

8  10 

It  was  found  that  the  AKR  could  saturate  at  10  -10  times 

the  initial  (noise)  levels.  This  is  adequate  to  produce  the  observed 

levels  of  AKR.  Almost  all  of  this  energy  comes  from  the  beam,  and  the 

4/3 

saturation  amplitude  was  seen  to  scale  as  ,  so  that  the  growth 

rate  of  the  beat  wave  determines  the  total  energy  the  electromagnetic 
wave  can  absorb.  Finally,  it  was  found  that  for  sufficiently  large 
density  cavities  with  growth  regions,  the  saturation  level  of  AKR  is 
relatively  independent  of  the  initial  level  of  the  EIC  density  fluctuations. 
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